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ABSTRACT 

The dispersion and mean trends of r-process abundances in metal-poor stars are discussed based on a 
model of diverse supernova sources for the r-process. This model is unique in that its key parameters are 
inferred from solar system data independent of stellar observations at low metallicities. It is shown that 
this model provides a good explanation for the observed dispersion and mean trend of Eu abundances 
over —3 <^ [Fe/H] <^ —1. It is also shown that this model provides a means to discuss r-abundances in 
general. For example, the Ag abundance in any metal-poor star with observed Eu and Fe abundances 
can be calculated from the model. This approach is demonstrated with success for two stars and can be 
further tested by future Ag data. The dispersion and mean trend of Ag abundances in metal-poor stars 
are also calculated for comparison with future observations. 

Subject headings: Galaxy: evolution — Stars: abundances — Stars: Population II 
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1. INTRODUCTION 

This paper discusses the dispersion and mean trends 
of r-process abundances (r-abundances) in "metal-poor" 
stars (with —3 ^ [Fe/H] <^ —1) based on a model of di- 
verse supernova sources for the r-process inferred from 
meteoritic data. Wasserburg, Busso, & Gallino (1996) 
showed that meteoritic data on the ^^^I and ^*^Hf in- 
ventory in the early solar system require two distinct fre- 
quencies for replenishment of the r-nuclei below and above 
mass number A ^ 130 in the interstellar medium (ISM). 
The frequency required for replenishment of the heavy 
r-nuclei {A ^ 130) is ~ (lO''' yr)"^ while that for the 
light r-nuclei {A <^ 130) is ~ 10 times less. These fre- 
quencies are consistent with the occurrence of Type II 
supernovae (SNII) in an average Galactic ISM (Wasser- 
burg et al. 1996; see also §2). The possible realization 
of different r-process production in two distinct kinds of 
SNII sources, the Ti (high-frequency, mainly producing the 
heavy r-nuclei) and C (low-frequency, mainly producing 
the light r-nuclei) events, was discussed by Qian, Vogel, & 
Wasserburg (1998). Based on the solar system data, Qian 
& Wasserburg (2000, QW) calculated the r-abundances 
resulting from a single Ti ot C event. Wasserburg & Qian 
(2000, WQ) showed that the results of QW were in general 
agreement with existing observations of metal-poor stars. 

The approach of QW and WQ is further developed here 
to study the dispersion in the abundance of a given r- 
element at a given metallicity, with the aim to gain some 
new insights into the evolution of r-abundances relative to 
Fe and to test SNII as the r-process site. The element Eu 
or Ag is chosen as the typical heavy or light r-element, re- 
spectively, for discussion. It is shown in §2 that the model 
of diverse SNII sources for the r-process inferred from me- 
teoritic data provides a good explanation for the observed 
dispersion and mean trend of Eu abundances in metal- 
poor stars. It is also shown in §3 that this model provides 
a means to discuss r-abundances in general. For example, 
the Ag abundance in any metal-poor star with observed 
Eu and Fe abundances can be calculated from the model. 



This approach is demonstrated with success for two stars 
and can be further tested by future Ag data. The dis- 
persion and mean trend of Ag abundances in metal-poor 
stars are also calculated for comparison with future ob- 
servations. Refinements of the model are discussed and 
conclusions given in §4. 

2. EU ABUNDANCES IN METAL-POOR STARS 

Consider the following simple picture for chemical en- 
richment of the Galactic ISM by SNII (e.g., Qian 2000). 
The ejecta of any individual SNII is mixed with a fixed 
amount of ISM corresponding to the total mass swept up 
by an SNII remnant, M^i^ « 3 x lO'' Mq (e.g., Thornton 
et al. 1998). The SNII rate per unit mass of gas f^^ /M^^ 
is assumed to be constant over Galactic history. Conse- 
quently, an average ISM in the Galaxy is enriched in the 
nucleosynthetic products from SNII at a frequency 
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where /q /Afgas is estimated using quantities for the 
present Galaxy. The frequency of Ti. events inferred from 
the meteoritic data is fn ~ /^bc ^^'^ that of £ events 
is fc ^ 0.1/^^. For definiteness, fn = (10^ yr)~^ and 
fc ~ (10^ yr)~^ are used below. It is further assumed 
that all SNII are either Ti, oi C events, with the average 
fraction of each kind being q — fn/{fn + fc) = 10/11 for 
Ti events or 1 — g = 1/11 for C events. The nucleosynthetic 
yields in an or £ event are taken to be constant. 

Over the Galactic history of « 10^° yr before solar sys- 
tem formation, w 10^ Ti and w 10^ C events contributed to 
the solar abundances. As the Ti events are mainly respon- 
sible for the heavy r-elements such as Eu, the Eu abun- 
dance resulting from a single Ti event must be (Eu/H)-^ ~ 
lO~^(Eu/H)0_j., where Eu/H is the number ratio of Eu to 
hydrogen and (Eu/H)0_r ~ 3.29 x 10"-^^ is the solar r- 
process value (Arlandini et al. 1999). In the spectroscopic 
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notation, logeT^(Eu) = log(Eu/H)7^ + 12 « -2.48. By ig- 
noring possible Eu contribution from the L events, the Eu 
abundance in a star can be considered as resulting from a 
number n-yi of H. events: 



loge(Eu) « loge7^(Eu) + logTi-^. 



(2) 



For example, HD 122563, HD 115444, and CS 22892-052 
with [Fe/H] = -2.74, -2.99, and -3.1 have loge(Eu) = 
-2.59, -1.63, and -0.93, respectively (Westin et al. 2000; 
Sneden et al. 2000), corresponding to nu « 1, 7, and 35. 

The three stars mentioned above all have [Fe/H] « —3 
but very different loge(Eu) values. In fact, this large 
dispersion in loge(Eu) at [Fe/H] « —3 has been demon- 
strated by more extensive observations (Mc William et al. 
1995; Burris et al. 2000). Application of equation (||) to 
the Eu data at [Fe/H] « —3 led WQ to conclude that 
very little Fe can be produced by the 7i events. It was 
shown that « 1/3 of the solar Fe inventory was provided 
by SNII (e.g., Timmes, Woosley, & Weaver 1995; Qian 
& Wasserburg 2001a). This Fe must be assigned to the 
« 10^ £ events that contributed to the solar abundances. 
Consequently, the Fe abundance resulting from a single 
£ event must be (Fe/H)^ « (Fe/H)0/3OO, corresponding 
to [Fe/H]£ = log(Fe/H)£-log(Fe/H)0 « -2.48. Further- 
more, WQ concluded that the Fe inventory at [Fe/H] w —3 
must be provided by the first very massive stars (with 
masses of J> IOOMq) before any SNII could have oc- 
curred. These stars may have produced other elements 
along with Fe, thus providing an "initial" inventory of met- 
als to the ISM. However, data on Ba (a heavy r-element) 
at -4 < [Fe/H] < -3 (McWiUiams et al. 1995; McWiUiam 
1998) show that the initial inventory of heavy r-elements 
is negligible compared with the contribution from a single 
H. event (WQ). Therefore, the initial Eu inventory can be 
ignored in equation (|) for loge(Eu). The initial inventory 
of primary concern here is that of Fe, hereafter denoted 
as [Fe/H]o — —3. As Type la supernovae (SNIa) con- 
tributed Fe only at times corresponding to [Fe/H] > — 1 
(e.g., Timmes et al. 1995; Qian & Wasserburg 2001a), the 
Fe abundance in a metal-poor star with —3 <^ [Fe/H] <J — 1 
can be considered as resulting from a number nn of L 
events: 



[Fe/H] w log (l0P<=/H]o + X 10P<=/"1^ 



(3) 



The range — 3 [Fe/H] < — 1 of interest here corresponds 
to nc ~ 0-30. 

Based on the above discussion, the Eu and Fe abun- 
dances in a metal-poor star can be considered as resulting 
from a mixture of 7i and C events. Given the average 
fractions of H, and L events among all SNII, the distri- 
bution of loge(Eu) expected for metal-poor stars with a 
given [Fe/H] can be calculated. Consider the occurrence 
of SNII in an ISM with the average fraction of Ji events 
being q. The probability to have a number n-yi of 7i events 
mixed with a number rin of L events is the same as that 
for the (n-H ^nc + l)th SNII to be the {nc + l)th L event. 
The latter is given by the so-called negative binomial dis- 
tribution: 



P{nn,nc) = 



{n-H + ncy. 



In particular, P{nT-i, n£ = 0) — (1 — q)q"'" . Thus the mix- 
tures of {nfi,nc) ~ (1,0), (7,0), and (35,0) corresponding 
to the Eu and Fe abundances in HD 122563, HD 115444, 
and CS 22892-052 would occur with probabihties of w 8%, 
5%, and 0.3% at [Fc/H] « -3 for q = 10/11. In general, 
among all the mixtures with a given nc the fraction having 
n-j-c < n^ is 



P{nH,nc). 



(5) 



nn=0 



Consequently, a fraction x of all the stars with [Fe/H] given 
by equation (m for the same nc as in equation (^ should 
have loge(Eu) < loge'^(Eu) « loge-H(Eu) -I- logn^. An 
"average" value of loge(Eu) at a given [Fe/H] can be de- 
fined as (loge(Eu)) w loge>^(Eu) -I- log(n-^), where 



1 



-inc + 1) 



(6) 



(l-g)"^+ig"«. 
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is the mean value of n-u for a given nc according to the 
distribution in equation (^. 

By taking logew(Eu) = -2.48, [Fe/H]o = -3, 
[Fc/H]£ = -2.48, and q = 10/11, the contours of 
loge^(Eu) are shown in Figure 1 as functions of [Fe/H] 
for X « 1%, 10%, 90%, and 99%. These can be inter- 
preted as the 99% (solid hues, x « 1% and 99%) or 90% 
(dashed lines, x w 10% and 90%) "confidence contours" 
for Eu abundance at a given [Fe/H]. If the model presented 
here is correct, then it should be very unlikely to find stars 
with values of loge(Eu) and [Fe/H] outside the 99% con- 
fidence contours. The dot-dashed line in Figure 1 shows 
(log e(Eu)) as a function of [Fe/H] and can be taken as the 
"mean" trend for evolution of Eu abundance relative to 
Fe. The Eu data of McWiUiam et al. (1995), Westin et al. 
(2000), Sneden et al. (2000), and Burris et al. (2000) are 
also shown in Figure 1. The data at —3 <^ [Fe/H] ^ — 1.3 
appear to be described quite well by the confidence con- 
tours and mean trend calculated from the model, espe- 
cially when typical observational errors of !=a ±0.1-0.2 for 
loge(Eu) are taken into account. As the model only con- 
siders Fe addition from the C events to the initial Fe inven- 
tory, it tends to overestimate nc [hence and log e^(Eu)] 
near the onset of SNIa Fe contribution at [Fe/H] ~ —1. 
This explains the few data points that are significantly 
below the contour for x « 1% at [Fe/H] ^ —1.3. These 
data points are expected to be inside the 99% confidence 
contours when SNIa Fe contribution is allowed for. The 
vertical segments of the contours for x « 1% and 10% are 
caused by rt^ = [loge^(Eu) = — cx)] for nc ~ I and 0, 
respectively. These segments will be modified by a refined 
model that takes into account a small initial Eu inventory 
associated with [Fe/H]o = — 3 and possible Eu contribu- 
tion from the C events (see §4). 

3. AG ABUNDANCES IN METAL-POOR STARS 

As discussed in §2, the model of diverse SNII sources 
for the r-process based on the characteristics of TC and C 
events inferred from meteoritic data provides a good expla- 
nation for the observed dispersion and mean trend of Eu 
abundances in metal-poor stars. There may be other possi- 
ble models that can achieve the same by separating Eu and 
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Fe production in a similar way. However, the model pre- 
sented here is unique in that its key parameters are inferred 
from solar system data independent of stellar observations 
at low metallicities. Furthermore, this model provides a 
means to discuss r-abundances in general. For example, 
the abundance of Ag (a light r-element) in a metal-poor 
star can be considered as resulting from a mixture of Ti 
and C events in addition to a possible initial Ag inventory 
associated with [Fe/H]o = —3. The numbers n-^ and nc 
for the mixture can be determined from the observed Eu 
and Fe abundances as discussed in §2. Given the initial Ag 
inventory logeo(Ag), the Ag abundance resulting from a 
single Ti. event loge7^(Ag), and that from a single C event 
loge£(Ag), the Ag abundance in a metal-poor star can be 
calculated as 



loge(Ag) sa log 



+ nc X ioi°s<^^(Ag) 



(7) 



Stars with [Fe/H] w —3 correspond to mixtures with 
n£ ?s (see eq. pi ) . If Ag abundances are known for two 
such stars with different observed Eu abundances (hence 
different n^, e.g., HD 115444 and CS 22892-052), then 
the parameters logeo(Ag) and loge-H(Ag) can be obtained 
from two equations similar to equation (^. However, so 
far the Ag abundance in only one star with [Fc/H] « —3 
has been published [loge(Ag) = -0.80 for CS 22892-052, 
Sneden et al. 2000]. For the discussion to proceed, it 
is assumed that the initial Ag inventory can be ignored 
[logeo(Ag) = -oo]. With {nn,nc) « (35, 0) for CS 22892- 
052 (see §2), loge7^(Ag) « —2.34 is obtained. As the solar 
r-abundance of Ag resulted from w 10^ Ti. and ~ 10^ C 
events (see §2), an equation similar to equation (|^) also ap- 
plies to the sun. The total solar photospheric abundance 
loge0(Ag) = 1.08 (Crawford et al. 1998) and the solar s- 
process abundance (Ag/H)© s = 3.45 x 10^^^ (Arlandini et 
al. 1999) correspond to logeQ,r(Ag) « 0.93, which gives 
loge£(Ag) « —1.40 (the Ag yield in an C event is w 9 
times higher than that in an H. event, as expected from 
the characteristics of these events). With the parameters 
logeo(Ag), loge-H(Ag), and loge£(Ag) chosen above, the 
Ag abundance in any metal-poor star with observed Eu 
and Fe abundances can be calculated by using equations 
§, (|), and (0). 

Crawford et al. (1998) reported Ag abundances or up- 
per hmits for 7 stars with [Fe/H] « -2.46 to -1.32, two 
of which also have observed Eu abundances from Burris 
et al. (2000). For these two stars the Ag abundances cal- 
culated from the above model are in good agreement with 
the observed values (see Table 1; the observational error 
for loge(Ag) is w ±0.2). The model can be further tested 
by future Ag data. As part of the test, stars that are pre- 
dicted to have high Ag abundances can be selected first for 
observation/analysis. For stars with observed Ag abun- 
dances but no Eu data. Figure 2 shows the dispersion and 
mean trend of Ag abundances over —3 < [Fe/H] ^ — 1 cal- 
culated as in the case of Eu (see §2) for logeo(Ag) = — oo, 
logeH(Ag) = -2.34, loge£(Ag) = -1.40, and q = 10/11. 
The Ag data of Crawford et al. (1998) and Sneden et 
al. (2000) are also shown in Figure 2. While the available 
data are fully consistent with the model, extensive Ag data 
for metal-poor stars are needed for a detailed comparison. 



The vertical segments of the confidence contours in Fig- 
ure 2 are caused by the assumption of logeo(Ag) = —oo 
and will be modified when the actual value of logeo(Ag) 
is known (see §4). 

4. DISCUSSION AND CONCLUSIONS 

The dispersion and mean trends of r-abundanccs in 
metal-poor stars have been discussed based on a model 
of diverse SNII sources for the r-process. This model is 
unique in that its key parameters are inferred from so- 
lar system data independent of stellar observations at low 
metallicities. It has been shown that this model provides 
a good explanation for the observed dispersion and mean 
trend of Eu abundances over — 3 ^ [F^/H] —1. The 
calculated confidence contours and mean trend for Eu in 
Figure 1 and those for Ag in Figure 2 represent a good 
test to distinguish SNII from neutron star mergers as the 
r-process site. Extremely high Eu and Ag abundances 
compared with the confidence contours in Figures 1 and 
2 would be expected for metal-poor stars if neutron star 
mergers were the r-process site (Qian 2000). If the results 
presented here are borne out by more extensive Eu and Ag 
data, then it is quite secure to conclude that SNII are the 
major sources for the r-process. Furthermore, the model 
presented here provides a means to discuss r-abundances 
in general. For example, the Ag abundance in any metal- 
poor star with observed Eu and Fe abundances can be cal- 
culated from the model. This approach has been demon- 
strated with success for two stars (Table 1) and can be 
further tested by future Ag data. The full approach for a 
large number of r-elements will be discussed in a separate 
paper (Qian & Wasserburg 2001b). 

As emphasized above, the parameters for the model pre- 
sented here are inferred from solar system data. In par- 
ticular, the fraction q of H events among all SNII is esti- 
mated from the frequencies of Ti. and C events suggested 
by the meteoritic data on ^^^I and ^*^Hf. While the spe- 
cific choice of 5 = 10/11 appears to describe the existing 
observations quite well [this is confirmed by the indepen- 
dent work of Fields, Truran, & Cowan (2001)], the actual 
value of q may be obtained from e.g., an extensive set of Eu 
data at [Fe/H] « —3. Stars with [Fe/H] f» —3 contain no 
contributions from the £ events. The probability for such 
stars to have received contributions from a number n-u of 
Ti. events is simply proportional to g"" (§2). Abundances 
of Eu in these stars can be converted into a histogram on 
the frequency of occurrences for different uh- Then the 
value of q can be obtained by an exponential fit to the 
histogram. It can be seen from Figure 1 that the exist- 
ing Eu data at [Fe/H] « —3 are rather inadequate for this 
purpose. Future observations in this region of [Fe/H] are 
highly desirable. 

The initial Eu inventory associated with [Fc/H]o = — 3 
and possible Eu contribution from the C events have been 
ignored. This gives rises to the vertical segments of the 
confidence contours in Figure 1. The actual initial Eu 
inventory may be estimated from the Eu abundances at 
[Fe/H] < —3 corresponding to the epoch for production 
by the first very massive stars. The actual Eu produc- 
tion in an C event is more complicated to determine. One 
possibility is to detect overabundances of r-elements on 
the surface of a binary companion to a neutron star (Qian 
2000). If the SNII producing the neutron star was an £ 
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event, then large overabundances of light r-elements such 
as Ag should be observed. Any associated overabundance 
of Eu, if detected, can be used to estimate the Eu produc- 
tion in an C, event. 

The initial Ag inventory have also been ignored by tak- 
ing logeo(Ag) = — oo. The actual value of logeo(Ag) 
can be determined when the Ag abundance in another 
star with [Fe/H] « -3 (e.g., HD 115444) is known. The 
upper limit on the Ag abundance in HD 140283 with 
[Fe/H] = -2.46 (Crawford et al. 1998) indicates that 
logeo(Ag) <j —1.15 (see Figure 2). In any case, extensive 
Ag data are needed to test and refine the model presented 



here. 
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Fig. 1. — The dispersion and mean trend of Eu abundances in metal-poor stars calculated from the model axe compared with data (open 
squares: McWilliam et al. 1995; asterisks: Westin et al. 2000; filled diamond: Sneden et al. 2000; filled circles: Burris et al. 2000). The solid 
and dashed lines correspond to the calculated 99% and 90% confidence contours for Eu abundance at a given [Fe/H]. The dot-dashed line is 
the calculated mean trend for evolution of Eu abundance relative to Fe. 

Fig. 2. — Same as Figure 1, but for Ag. Open triangles (upper limits) and filled squares are new symbols indicating data from Crawford et 
al. (1998). 



Table 1 

Comparison of calculated Ag abundances with data 



Star 


[Fe/H] 


nc 


log e(Eu) 




logeobs(Ag) 


logecal(Ag) 


HD 2665 


-1.95 


3 


-1.23 


18 


-0.88 


-0.70 


HD 6755 


-1.60 


7 


-0.56 


83 


-0.16 


-0.18 



Note. — Data on Fe and Ag [log eobs (Ag)] are taJsen from Crawford et al. 
(1998), and data on Eu from Burris et al. (2000). The log ecal (Ag) values are 
calculated from the model (see text). 



